Spin freezing transition and non-Fermi-liquid self-energy in a 3-orbital model 
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A single-site dynamical mean field study of a three band model with the rotationally invariant 
interactions appropriate to the t2 g levels of a transition metal oxide reveals a quantum phase transi- 
tion between a paramagnetic metallic phase and an incoherent metallic phase with frozen moments. 
The Mott transitions occurring at electron densities n = 2, 3 per site take place inside the frozen 
moment phase. The critical line separating the two phases is characterized by a self energy with the 
frequency dependence E(w) ~ yfu and a broad quantum critical regime. The findings are discussed 
in the context of the power law observed in the optical conductivity of SrRuOs. 
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The 'Mott' metal-insulator transition plays a central 
role in the modern conception of strongly correlated ma- 
terials [H, dj . Much of our understanding of this transi- 
tion comes from studies of the one-band Hubbard model. 
Here, the transition is generically masked by antiferro- 
magnetism, but if this is suppressed (physically, by in- 
troducing lattice frustration or mathematically, by ex- 
amining an appropriately restricted class of theories such 
as the paramagnetic-phase single site dynamical mean 
field approximation 3|) a transition from a paramagnetic 
metal to a paramagnetic insulator is revealed. The prop- 
erties of the paramagnetic metal phase near the tran- 
sition play a central role in our understanding of the 
physics of correlated electron compounds. 

While one band models are relevant to many materi- 
als including the high temperature superconductors and 
some organic compounds, many systems of interest in- 
volve multiple correlated orbitals for which the physics 
is richer and less fully understood. Multiorbital models 
have been studied in Refs. @, H, % 0, B [3, [13]. New 
physics related to the appearance of magnetic moments 
has been considered in the context of the orbitally se- 
lective Mott transition which may occur if the orbital 
degeneracy is lifted [!l G3, El El [lit , but for orbitally 
degenerate models it seems accepted that the essential 
concepts of a paramagnetic metal to paramagnetic insu- 
lator transition and a strongly correlated paramagnetic 
metal phase can be carried over from studies of the one- 
band situation. 

In this paper we show that this assumption is not cor- 
rect. We use the single-site dynamical mean field approx- 
imation to demonstrate the existence of a quantum phase 
transition between a paramagnetic Fermi liquid and an 
incoherent metallic phase characterized by frozen local 
moments (a spin-spin correlation function which does not 
decay to zero at long times). We show that for densities 
per site n — 2,3 the Mott transition occurs within or 
at the boundary of the frozen moment phase. As Costi 
and Liebsch have noted in the context of an orbitally 
selective Mott system, the presence of frozen moments 



may be expected to influence the Mott transition [15| . 
The new phase appears for multiple orbitals, a differ- 
ent number of electrons than orbitals and a rotationally 
invariant on-site exchange J7/3 > J > 0. While this sit- 
uation is relevant to many compounds, it has has only 
recently become amenable to systematic study follo wing 
the development of improved exact diagonalization [161 ] 
and efficient Monte Carlo methods 17], [l8| ■ 

The transition into the frozen moment phase produces 
a non-Fermi liquid self energy which varies as y/cJ, and 
the quantum critical regime of the phase transition is suf- 
ficiently large that the non-Fermi liquid behavior is ob- 
servable over a wide range of temperature, interaction 
strength, spin polarization and carrier concentration. 
The non- Fermi- liquid self energy we find bears an intrigu- 
ing resemblance to the self energy inferred from optical 
conductivity measurements on SrRu0 3 [3 S3, El . 

We study a model inspired by titanate, vanadate and 
ruthenate transition metal oxide compounds, in which 
the relevant electrons reside in ^g-symmetry d-orbitals. 
In the commonly occurring pseudocubic structure the 
t2 g levels are three-fold degenerate, and the formal va- 
lence is such that the levels contain 1 electron (Ti com- 
pounds), 2 electrons (some V compounds) or 4 electrons 
/ 2 holes (Ru compounds). The same physics occurs 
for the two-orbital models relevant to the nickelates and 
other e s -based systems. The electronic Hamiltonian is 

H = ffband + Si-^loc) witn #band = E«,p £ p ^p,t^p,<r 

and the local Hamiltonian on each site i given by 



(1) 



+ 51 U'n a , a np,- a + (U' - J)n a . a np. a 

a>j3,a 
a=£/3 

Here, a = 1, 2, 3 is the orbital index, a =|, J, the spin in- 
dex, fx is the chemical potential, U (U') the intra-orbital 
(inter-orbital) Coulomb interaction, and J the coefficient 
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FIG. 1: Phase diagram for J/U = 1/6 and fit = 50, 100 in the 
space of density n and interaction strength U. The light line 
with circles or diamonds indicates a phase transition between 
a Fermi liquid metal and a "frozen-moment" metal. The black 
lines mark the regions of Mott insulating behavior. 



of the Hund coupling. We adopt the conventional choice 
of parameters, U' = U — 2 J. 

To study the model we use the single-site dynami- 
cal mean field approximation 0] which ignores the mo- 
mentum dependence of the self-energy and reduces the 
original lattice problem to the self-consistent solution of 
a quantum impurity model given by the Hamiltonian 
Hqi = H\ oc + H hyb + i?bath with H hyh describing the 
exchange of electrons between impurity and bath. Our 
data were computed for a semi-circular density of states 
with band-width At (so the model is particle-hole sym- 
metric about the density n = 3) , using the continu- 
ous time QMC solver of Ref. [T3, EH]. We investigate 
the electron self energy and the imaginary-time 

impurity-model spin-spin and orbital-orbital correlators 
(O(r)O(O)) with O representing either the electron spin 

density S z = lJ2 a |( rf l,T rf a>T ~ d lt,l da >^ or tne orbital 
density h a — <P a a d a , a - Attention is restricted to 
solutions which preserve spin and orbital rotational sym- 
metry at the level of the impurity model. 

Figure [T] presents our calculated phase diagram in the 
space of density n and interaction strength U for the 
ratio J/U = 1/6. The Mott insulating phases of the 
model are shown as heavy solid lines. The light line with 
circles or diamonds is our new result: a phase boundary 
separating a small n small U Fermi liquid phase from a 
frozen moment phase at larger n and larger U. Other 
values of < J/U < 1/3 give similar results. For J = 
the new phase does not exist while for J > U/3 the term 
U 1 — J = U — 3J changes sign and the physics of the 
model becomes different. 

We may define the phase boundary using the im- 
purity model spin-spin correlation function Css{ T ) — 
(S z (t)S z (0)} 7 shown in the upper panel of Fig. [2] for 
U/t = 8 and several values of n. In a Fermi liquid at low 
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FIG. 2: Upper panel: imaginary time dependence of the 
spin-spin correlation function (S z (0)S z (t)) (positive correla- 
tion function, full symbols) and orbital correlation function 
(m (0)n2(i~)} (negative correlation function, open symbols) for 
U = St and carrier concentrations n indicated. Lower panel: 
variation with doping of the temperature dependence of the 
spin-spin correlation at r = 13/2. The error bars are large at 
smaller n because the midpoint spin correlator is very small. 
The black line indicates the n-value of the phase transition 
deduced from an analysis of the self energy. 



temperature T, Css(t) ~ (T/ sin(7rrT)) 2 for imaginary 
times t sufficiently far from either r = or t = l/T. Our 
results are consistent with this form in the Fermi liquid 
phase, but in the non- Fermi- liquid phase Css is seen to 
approach a constant at long times indicating the presence 
of frozen moments. We also plot in Fig.[^]the correspond- 
ing orbital correlation function, which is seen to decay 
rapidly with time on both sides of the phase transition. 
For a more quantitative analysis we studied the temper- 
ature dependence of Ci/ 2 = Css{ T = ^f)- In a Fermi 
liquid, C1/2 ~ T 2 while in the frozen moment phase C1/2 
becomes temperature independent at sufficiently low T. 
Within our numerical accuracy, we find that at the transi- 
tion point, Css(t = 27) ~ T. The lower panel of Fig. [2] 
shows how the ratio C 1/2 (T = 0.02t)/C 1/2 (T = OMt) 
changes from the value 4 expected in the Fermi liquid 
phase to the value 1 expected in the frozen moment 
phase. 
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FIG. 3: Doping dependence of the imaginary part of the Mat- 
subara axis self energy for U/t = 8, /3t = 50 (full symbols) 
and 100 (open symbols) at dopings indicated. The dashed 
lines are proportional to (oj n /t) a , a — 0.5, 0.54, 0.62 (from 
top to bottom). 
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FIG. 4: Results of a fit of the computed Matsubara-axis 
self energy to the scaling form — ImE/t — C + A(uj n /t) a at 
temperatures T = t/50 (diamonds; blue on-line) and T = 
t/lQQ (circles, red on-line) at density n — 2 and interaction 
strengths indicated. The plot illustrates the wide quantum 
critical regime of the spin freezing transition. 



The phase transition has consequences for the electron 
self energy E(w) . In a Fermi liquid at low T the imaginary 
part of the real axis scattering rate S (u) ~ max(w 2 , T 2 ) 
so the imaginary part of the Matsubara axis self energy 
Yi(iui n ) ~ iuj n at small u> n . Frozen moments may be 
expected to scatter electrons so that the real axis self 
energy is £ (lu) ~ T, implying on the Matsubara axis 
lnaS(iu) n — > 0) = iTsgn(uj n ). At the critical point we 
expect a power law behavior £ (u>) ~ lu 01 ; if a < 1 then 
S(«o; n — > 0) ~ (iu} n ) a . Figure [3] shows that the imag- 
inary part of our computed Matsubara-axis self energy 
is consistent with this behavior, vanishing linearly in the 
Fermi liquid phase, exhibiting an intercept in the frozen 
moment phase and an exponent a ~ 1/2 at n c . The be- 
havior of the self energy in the region not too far from 
the transition is described by a quantum critical crossover 
function. Our resolution is not sufficient to identify this 
function precisely. We have fit the self energy to the 
form lmE(u> n )/t = C + A(ui n /t) a , recognizing that the 
varying value of a represents the different regimes of the 
crossover function. In Fig. [4] we plot the exponent a and 
intercept C extracted from the fit as the transition line 
is crossed by varying interaction strength at fixed den- 
sity n = 2. Plots of this sort enable a simpler and more 
accurate location of the transition line than an analysis 
of Cssij) and were used to compute the phase boundary 
shown in Fig. [TJ 

Figures [3] and [¥] indicate that an approximately square 
root behavior of the self energy persists for a range of 
carrier concentrations and interaction strengths near the 
critical line. To further investigate the effect of frozen 
spins on S we have studied the magnetic field depen- 
dence of the self energy. As the spins are polarized by 
the field, the intercept decreases and the apparent power 
law changes towards the Fermi liquid exponent of a = 1 . 
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FIG. 5: Self energy for the majority spin as a function of 
magnetic field at U/t — 7, fit = 50 and n = 2, slightly above 
the critical point for the glass transition. 



However, near the critical point (see Fig. [5]) a changes 
slowly and the approximate square root frequency depen- 
dence remains visible even for relative magnetizations of 
the order of 50%. 

We speculate that the ^fw self energy found here is 
related to the still-mysterious optical conductivity er(fi) 
of SrRu0 3 and CaRu0 3 0, Sfl [El- These materials 
are pseudocubic perovskites with two holes in the t2 9 d 
shell and their behavior should be roughly modelled by 
our 3-orbital calculation at carrier density n = 2 (we 
neglect effects such as the buckling of octahedra or spin- 
orbit coupling). In these materials a(Q) varies approxi- 
mately as l/^/fi at higher frequencies, but rolls over (in 
a way which depends on temperature and on the ma- 
terial) to an approximately constant behavior. In the 
momentum-independent self energy approximation, ver- 
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tex corrections may be neglected and if the self energy is 
small compared to the bare electron bandwidth the real 
part of the optical conductivity may be written (/ is the 
Fermi function) 



cr(ft) 



dio f(uj + ft) - f(uj) 

~n- s(w + ft) + 



(2) 



so that a square-root self energy implies that the conduc- 
tivity varies as 1/ \/ft for frequencies such that S(ft) > ft. 
Our finding, that the phase boundary runs close to n = 2 
for a range of interaction strengths suggests that an ex- 
treme degree of fine tuning is not necessary to produce 
a square root self energy. It is important to note that 
in SrRuC>3 the square root behavior persists well into 
the ferromagnetically ordered phase at low temperature 
whereas in CaRu03 the conductivity develops a non- 
monotonic structure at low ft and higher T . 

The spin-freezing transition itself is of theoretical in- 
terest. Studies, to be presented elsewhere, of the model 
without the DMFT self consistency condition show that 
the transition exists at the level of the impurity model. 
We have also studied the degree to which the different 
eigenstates of H\ oc are represented in the partition func- 
tion. We find that for J > at couplings (U > 4t) only a 
few states are relevant. The large [/-density driven tran- 
sition is marked by a change in the dominant states from 
the one-electron states S = 1/2, L = 1 to a, nine- fold de- 
generate manifold of two electron states with S = 1 and 
L = 1, with the two manifolds becoming degenerate at 
the transition. The interaction-driven transition is on the 
other hand marked by a change in the weight of the two 
subleading states S = 1/2, L = 1 and S = 3/2 L = 0, im- 
plying a change in the magnitudes of coupling strengths. 
The orbital degree of freedom is important to stabilize 
the metallic phase at relevant interaction strengths (the 
two orbital model with two electrons and J/U = 1/6 
is insulating for U > 3.7i [3])- Trying to suppress the 
L = 1 orbital angular momentum states by applying a 
crystal field rapidly leads to an insulator. 

In conclusion, we have shown that in a model, rele- 
vant to transition metal oxides with partly filled d-shells, 
with several electrons in a three-fold degenerate level, an 
apparent spin freezing transition occurs. While it is pos- 
sible that the effects could be due to a rapid decrease of 
the spin coherence scale to values below the range ac- 
cessible to us, the square root self energy and T-linear 
spin-spin correlation function are strong evidence for an 
actual T = transition. The frozen moment phase re- 
sults from a calculation in which spin rotation symmetry 
was enforced and may thus be preempted by a state with 
broken spin rotation symmetry. The exact nature of the 
broken symmetry state will depend upon the underlying 
lattice (for the models studied here neither ferro nor an- 
tiferromagnetic solutions are stable; we expect that the 
true ground state is either an incommensurate magnet 
or is phase separated). However, we emphasize that in 



contrast to the situation in one-band models, a transi- 
tion remains even if long ranged order is suppressed in 
the calculation. 

Our results open many new research directions. The 
theory of the multielectron Mott transition must be re- 
examined, as our results suggest that even for models 
with high degrees of spin and orbital symmetry a para- 
magnetic metal-paramagnetic insulator transition gener- 
ically does not exist (at least for more than one elec- 
tron per site). The magnetic phase diagram of multior- 
bital models must be explored. We have shown that this 
new transition is driven by an apparently new class of 
impurity-model phase transitions which should be clas- 
sified and explored. The interplay with lattice-driven 
crystal field splittings and the orbitally selective Mott 
transition [3, EH is an open and interesting question. 
More generally the implications of the rich structure of 
the Slater-Kanamori interactions for the properties of the 
strongly correlated metal phase deserve further investi- 
gation. 

The calculations have been performed on the Hreidar 
beowulf cluster at ETH Zurich, using the ALPS-library 
0. PW and AJM acknowledge support from DOE-ER- 
16469 and EG and MT from the Swiss National Science 
foundation. 
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